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Long-lived Neel states in antiferromagnetic quantum spin chains with strong uniaxial 
anisotropy for atomic-scale antiferromagnetic spintronics 
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It has been experimentally established that magnetic adatoms on surfaces can be arranged to form 
antiferromagnetic quantum spin chains with strong uniaxial anisotropy and Neel states in such spin 
systems can be used to realize information storage. Here, we investigate eigen states, quantum spin 
dynamics, and life times of Neel states in short antiferromagnetic quantum spin chains with strong 
uniaxial anisotropy on the basis of numerical exact diagonalization method. We show rigorously 
that as long as the uniaxial anisotropy is very strong, the ground state and the first excitation state, 
being nearly degenerate, are safely separated from the other states and thus dominate the quantum 
dynamics of the Neel states. Through further numerical analysis, we achieve a powerful life-time 
expression of the Neel states for arbitrary spin and model parameters. It is interesting that for the 
famous Fe adatom chains on CU 2 N surface, 14 or 16 Fe adatoms are enough to obtain a practical 
long life-time for Neel state storage of information. These should be applicable to other similar 
antiferromagnetic spin systems for atomic-scale antiferromagnetic spintronics. 


1. Introduction 

It is inspiring that adatom-based antiferromagnets 
have been realized on semiconductor surfaces and used 
for novel magnetic information storage because their Neel 
states can be stabilized by strong uniaxial single-ion mag¬ 
netic anisotropyP-Q. Such nanomagnets can be fabri¬ 
cated adatom by adatom, and their spin anisotropy can 
be controlled P-ITJ. Thus, one can make antiferromag¬ 
netic chains, bi-chains, nano-ribons, or nano-sheets con¬ 
sisting of several or tens of adatom spins with strong 
magnetic anisotropy and adjustable inter-spin interac¬ 
tions. Spin chains are of much interest because they be¬ 
long to an important category of Heisenberg spin mod¬ 
els. In fact, various one-dimensional antiferromagnetic 
Heisenberg models have been intensively investigatedP- 
[Tsjl . For S' = 1, there exists an inetresting Haldane 
topological phase if there is no strong uniaxial magnetic 
anisotropy |la - [ld |. On experimental side, one usually 
use high spins with strong uniaxial single-ion magnetic 
anisotropy in adadtom spin systems [ll-l7|. It is known 
that strong uniaxial single-ion anisotropy is necessary 
to achieve stable Neel states. Experimentally, electrons 
currents injected through STM tips have been used to 
control the Neel states for information storage^. On 
theoretical side, some efforts have been made to under¬ 
stand and explore controlling the adatom-spin antifer¬ 
romagnets with spin-polarized electron current [l^ [^ . 
spin current . and mechanical oscillator (2^ and to in¬ 
vestigate symmetry effects on spin switching of single 
adatoms [^. It is believed that more significant ad¬ 
vances and deeper insight in this field can likely lead to 
an atomic-scale antiferromagnetic spintronics. 

Here, we investigate the intrinsic quantum dynamics 
and life times of Neel states in the quantum Heisen- 
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berg antiferromagnetic chain model consisting of 2N 
spins (S' > 1) with strong uniaixal single-ion anisotropy. 
We accurately calculate eigenvalues and eigenfunctions 
through exact diagonalization, and rigorously show that 
the ground state and the first excitation can be both 
safely separated from the other states and well de¬ 
scribed with the two Neel states as long as the single¬ 
ion anisotropy is very strong. Then, we thereby calcu¬ 
late the switching rates and life times of the Neel states. 
Surprisingly, we achieve a unified powerful expression of 
the life times through fitting our accurate numerical re¬ 
sults. More importantly, for the Fe-adatom spin antifer¬ 
romagnets on CU 2 N semiconductor surface[ll-Q, 2N= 14 
or 16 is large enough to achieve practical life times of Neel 
states for information storage. More detailed results will 
be presented in the following. 

2. Results and discussion 

2.1 Spin Model and eigenstates 

We start with general one-dimensional quantum 
Heisenberg antiferromagnetic model with strong uniax¬ 
ial single-ion anisotropy, 

2N-1 ^ ^ 2N 

H = jY.S.-S.+^-DY,{St)\ ( 1 ) 

i=l i=l 

where the total number of the spins is 2N, the parame¬ 
ter J (> 0) is the antiferromagnetic exchange constant, 
(> 0) is used to characterize the single-ion magnetic 

anisotropy in the z axis, and Si is the spin operator at 
site i, satisfying open boundary condition. Here, we do 
not need any inhomogeneous effective magnetic field to 
split the two Neel statesjl, [l^, but to experimentally 
prepare a specific Neel state, one can use an STM tip 
to inject a ^in-polarized electron current on the first 
adatom spin(3| . If being applied to similar spin rings with 
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FIG. 1. The eigen energies {E/J) of Hamiltonian ([T]) for S = 1 (a) and S = 3/2 (b), with D/J = 10, 3, and 1. 


TABLE I. The A^-dependent energy gap /S.E between the ground state G and the first excitation Ei for two spin values (1 and 
3/2) and three D/J (1, 3, and 10). 




A = 1 

to 

CO 

II 

A = 4 

A = 5 

S = 1 

D/J = 1 
D/J = 3 
D/J = 10 

0.561 

0.275 

9.48 X 10"^ 

0.158 

3.66 X 10"^ 
4.41 X 10“® 

4.76 X 10"^ 
4.95 X 10"® 
2.05 X 10"'‘ 

1.43 X 10"^ 

6.68 X 10"'‘ 

9.56 X 10"® 

4.32 X 10"® 
9.03 X 10"® 
4.46 X 10"^ 

S = 3/2 

D/J = 1 
D/J = 3 
D/J = 10 

0.228 

4.53 X 10"^ 
5.09 X 10"® 

1.24 X 10"^ 
4.44 X 10“^ 
5.69 X 10"® 

6.64 X 10"'‘ 
4.34 X 10"® 
6.36 X 10"® 

3.56 X 10"® 
4.25 X 10"® 
7.96 X 10"^^ 

1.90 X 10"® 
4.15 X 10"^° 
8.81 X 10"^® 


uniaxial anisotropy, such as antiferromagnetic molecule 
wheels [i^, the Hamiltonian (1) needs some modification 
to make the spin operators satisfy periodic boundary con¬ 
dition. 

Using Sf = Sf ± we have S, ■ S^+i = Sf ■ -h 
-I- S~ ■ The ideal Neel states (|iVi) and 

\N 2 )) are certainly not the eigenstates of Hamiltonian ([1]) 
due to the transverse part including the raising and low¬ 
ing operators 5/^, but the strong single-ion anisotropy 
D in the z axis makes the spin tend to orient in the 
z axis. Consequently, there are large {D/J dependent) 
weight of the two Neel states in the ground state and low 
excitations, which implies that the two Neel states for 
large D/J can be stable enough to be used for informa¬ 
tion storage^. Using exact diagonalization method}^ 
to the Hamiltonian (1), we can obtain the spin eigenval¬ 
ues and eigenfunctions. For convenience, we shall use J 
as our unit in the following, which means that anisotropy 


parameter D and energy E can be scaled in terms of J. 

In Fig. I we present the energy eigenvalues depending 
on N (1 through 5) and D/J (10, 3, and 1) for S' = 1 
and 3/2. For each of the cases, the ground state G and 
the first excitation Ei are both separated from the other 
states. The trend is that the separation increases with 
D/J and S. The corresponding energy gaps between 
the ground states and the first excitation ones are sum¬ 
marized in Table HI It is clear in the table that the gap 
decreases with A, S, and D/J. The weights of the Neel 
states in the ground state (G) and the first excitation 
(El) as functions of N are presented in Table HIl for the 
two spin values and the three D/J ones. It can be seen 
that the Neel weights increase with D/J, but decrease 
with N. Except the special case of A = 1, the Neel 
weights increase with S, too. Because we are interested 
in the cases with strong uniaxial anisotropy, the ground 
state and low excitations are far from the regime of the 
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Haldane state (T6l - [l8l| . 


TABLE II. The A^-dependent Neel weights in the ground state 
(G) and the first excitation (Ei) for different S and D/J. 





N = 1 

N = 2 

A = 3 

A = 4 

A = 5 

5=1 

D/J = 

= 1 

G 

0.864 

0.750 

0.693 

0.634 

0.575 




El 

1.000 

0.854 

0.735 

0.650 

0.580 


D/J = 

= 3 

G 

0.964 

0.935 

0.916 

0.892 

0.867 




El 

1.000 

0.958 

0.921 

0.893 

0.868 


D/J = 

10 

G 

0.996 

0.991 

0.988 

0.985 

0.980 




El 

1.000 

0.994 

0.989 

0.985 

0.980 

S = 3/2 

D/J = 

= 1 

G 

0.900 

0.875 

0.822 

0.770 

0.720 




El 

0.968 

0.882 

0.823 

0.770 

0.720 


D/J = 

= 3 

G 

0.985 

0.946 

0.957 

0.941 

0.926 




El 

0.991 

0.973 

0.957 

0.941 

0.926 


D/J = 

10 

G 

0.999 

0.995 

0.994 

0.992 

0.991 




El 

0.999 

0.996 

0.994 

0.992 

0.991 


It is easy to prove that the total spin z-component 
= J2i is conserved because it is commutable with 
the Hamiltonian (1). All the energy eigenstates can be 
classified in terms of the eigenvalue of S^. Generally 
speaking, for a finite antiferromagnetic chain with 2N 
spins, the ground state is a spin single state. When D/J 
is very large, the ground state can be approximately con¬ 
structed with a superposition of the two Neel states. For 
general D/J, we can always construct the following two 
eigenstates from the Neel states. 

|+) = ci(|Ai) + |iV 2 )+Oi|(!I+) + .-.) 

\-)=C2{\Ni)-\N2)+02\0.) + ■■■), 

where \0±) is defined as ‘ “^i+i + ^7 ' 

>?i+i)(|-^i)±|N 2 )), and Ci, C 2 , Oi, and O 2 are coefficients 
to be determined. Actually, our exact diagonalization re¬ 
sults show that when 2NS is even, the ground state G is 
|-|-) and the first excitation Ei is |—); and when 2NS is 
odd, we have G=|—) and Ei = |-|-). This is in accordance 
with the theoretical results obtained by spin coherent 
state path integralfl^ [l^- The higher excitation states 
with Ei {i > 2) can be constructed in the similar way. 

2.2 Quantum dynamics of Neel states 

We shall mainly focus on the subspace of the states 
with = 0 because the ground state and the low exci¬ 
tation states including the Neel states belong to this sub¬ 
space, and however, we shall turn to other states when 
we discuss the effect of temperature. Eor convenience, we 
shall use \g) and |ei) {i > 1) to denote all the eigenstates 
in the = 0 subspace. Because this subspace is closed 
under the Hamiltonian (1), the time evolution of the two 
Neel states can be expanded as 

Wait)) = + E (3) 

j>i 


where Ej and /“ (j > 0, a =1,2) are the eigenvalues 
and expansion coefficients of the j-th eigenstates. Here, 
of course, we have |iVa(0)) = \Na), | 5 )=G, and |ei)=Ei. 
Then, the weight of \Na) in \Na{t)) can be expressed as 

X^W = l/o (4) 

t>i 

where AEj = Ej —Eq- The total Neel weight of |Ai) and 
IA 2 ) in |iVi(t)) can be defined as WN{t) = Xi(t) 

WN(t) reflects how well the Neel states describe the quan¬ 
tum antiferromagnetic chain. The two-state approxima¬ 
tion results in a simplified expansion of \Na{t)), such as 

|Ai(t)) « cos(^t)IAi) + sin(^t)|Ar2). (5) 

We present xf(t) and W^it) in Eig. 2 for D/J=10 , 3, 
and 1. Eor the two-state approximation is also 

presented for comparison. It is clear that xf(t) is a pe¬ 
riodic function of t and W^it) is almost a constant ex¬ 
cept a narrowly oscillating noise due to the higher states. 
Eor small D/J such as 1, the maximal value of Xi(0 
is approximately 0.8 and the Neel weight W^it) is less 
than 0.9, but for large D/J such as 10, Xi(^) can be well 
described with cos^(t/2T) and the Neel weight becomes 
larger than 0.99. Here, the time period is equivalent to 
P = 2ttT, and 1/T reflects the switching rate (or fre¬ 
quency) between the two Neel states. It is surprising 
that for this case of S' = 2 and N = 2, T increases by five 
orders of magnitudes when D/J changes from 1 to 10. 

2.3 Life times of Neel states 

Because Xi(t) is a well-defined periodic function of 
t, the quantity T, the time spent by a switching cir¬ 
cle between the two Neel states, can be used to char¬ 
acterize the life times of the Neel states. In the case of 
two-spin chains {N = 1) with S < 3, we can calculate 
eigenstates and \Na{t)) exactly. For S = 1, we obtain 
AEi = jy4:{D/J)'^ +4:D/J + 9-2D-l]/2, and T can 
be expressed as {2D/J -I- l)/2 when D/J is large. For 
higher spins, we can achieve T oc {2D/J+1)‘^^~^ for both 
integer and half odd integer spins by using a usual per¬ 
turbation method. Generally speaking, we can also use 
exact diagonalization method to calculate AE and T for 
arbitrary S and N. In Fig. 3 we present the calculated T 
as functions of H/J for A = 1, with S taking nine values 
from 1 to 5. In Fig. 4 we present our accurate calculated 
T curves for N=l through 4 and 5'=1, 3/2, 2, and 5/2. 

It is very interesting that all these {D/J)-T curves can 
be satisfactorily fitted with one simple function, 

T = A^(2^ + 1)^(2^-L (6) 

where A is a constant depending on S and N only. It 
is surprising that, as we show in Fig. 5, A can be well 
fitted with A = ba^, and furthermore the parameters a 
and b can be well fitted with a = 0.2427 x 4.1545'^ and 
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FIG. 2. The time dependence of Xi(f) (black lines in the left 
column) and Wjv(t) (the right column) for 5 = 2, 2N — 4, 
and D/J = 10 (a,d), 3 (b,e), and 1 (c,f), respectively. The 
two-state approximated results of Xi (f) (blue dash lines in the 
left column ) are also presented for comparison. 


(a) (b) 



FIG. 4. The life times T (in unit of h/ J) of the 2N chains as 
function of D/J for 5 = 1 (a), 5 = 3/2 (b), 5 = 2 (c), and 
(d) 5 = 5/2, fitted with Eq. (5). 


b = 0.5007 X The fitted data of A, a, and b 

are summarized in Table IIIII Consequently, we obtain a 
unified expression for T as functions of D, J, N, and S. 
It can be used to extrapolate T with given D and J for 
higher S and larger N. It should be pointed out that 
although T increases with increasing D or decreasing J, 
too small J will be harmful to stability against thermal 
fluctuations. Although T increases exponentially with N 
increasing, one cannot use too long spin chains for prac¬ 
tical information storage because the Neel weight will 
decrease with N increasing. Therefore, for a practically 
useful system, one should keep a balance between a large 
T and a good stability of the Neel states. 



FIG. 3. The life times T (in unit of ft/J) of the two-spin 
chains as functions of D/J for nine 5 values, fitted with Eq. 
(5). 


TABLE III. Fitted results of A, a, and h for the antiferromag¬ 
netic spin-5 chains including 2N spins. 


Spin 

A 

a b 

N=1 N=2 N=3 A=4 

5=1 

0.506 0.519 0.531 0.542 

1.023 0.4952 

5=3/2 

0.445 0.909 1.839 3.725 

2.029 0.2202 

5=2 

0.500 2.039 8.527 36.47 

4.178 0.1183 

5=5/2 

0.640 5.245 47.49 371.9 

8.415 0.07596 

3=2, 

0.894 15.74 289.2 4871 

17.68 0.05080 


2.4 Long life times in real adatom spin chains 

For the short antiferromagnetic chains of Fe adatom 
spins on CU 2 N surface, experimental result reveals that 
S = 2, J = 0.7 meV, and D = 1.87 meV In this case, 
we have D/J = 2.67, belonging to the regime of strong 
uniaxial anisotropy. As a result, we obtain a simple for¬ 
mula of T (in second) depending on N, 

T = 1069^ X 1.120 X 10"^^. 


(7) 
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FIG. 5. (a) Fitting of the numeric results of A depending on 
N for different S'; (b) fitting of a depending on S; (c) fitting of 
b depending on S; and (d) Plot of the life time T (in second) 
depending on N defined in Eq. (6) for the Fe adatom spin 
chains on CU 2 N surface. 

We plot it in Fig. 5(d). The expression (6) implies that 
the switching rate (1/T) will be decreased approximately 
by 1000 times when we add two more Fe-adatom spins to 
the chain, which is consistent with the low-temperature 
limit of the experimental results Q. This implies that 
the life times T can be very long, reaching 1.9 days, 5.7 
years, and 6057 years when N is equivalent to 6, 7, and 
8, respectively. These results show that for such anti¬ 
ferromagnetic chains, 14 or 16 spins (for 2N) should be 
enough to achieve stable Neel states at ultra-low temper¬ 
atures (a few kelvin) for practical information storage. 


On the other hand, for practical usage, we need to 
consider other factors affecting the life times of Neel 
states. First, we consider possible transverse single-ion 
anisotropy -E and transverse magnetic field which ap¬ 
pear as additional — (5f)^] in the 

Hamiltonian. Our calculations reveal that as long as D/ J 
is not less than 1, there is little change in T even when 
Ej J and ^B^/J reach to 0.2. Then, we investigate effect 
of spin exchange anisotropy on T, showing that the effect 
is very small for D/J > 1. Therefore, our T results are 
robust and technically sound. 

3. Conclusions 

In summary, we have investigated the intrinsic quan¬ 
tum dynamics and life times of Neel states in the quan¬ 
tum Heisenberg antiferromagnetic chains with strong 
uniaixal single-ion anisotropy. For typical values of spin, 
chain length, and magnetic anisotropy, we have used ex¬ 
act diagonalization method to accurately calculate eigen¬ 
values and eigenfunctions, and shown rigorously that the 
ground state and the first excitation are both safely sep¬ 
arated from the other states and can be well described 
with the two Neel states as long as both D/J and S are 
large enough and N is not too large. Through investigat¬ 
ing accurate time evolution of the Neel states, we have 
determined their switching rates and life times. Surpris¬ 
ingly, we have achieved a unified powerful expression of 
the life times for arbitrary values of N, S, D, and J. 
Furthermore, we show that for the Fe-adatom spin anti- 
ferromagnets on CU 2 N semiconductor surface[lll3, 2N= 
14 or 16 is large enough to achieve practically long life 
times for the Neel states. These theoretical results should 
be useful to help realize the Neel state storage of infor¬ 
mation and atomic-scale antiferromagnetic spintronics. 
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